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唱 体 塑性 理论 对 单 唱 和 双 唱 的 压 痕 行 为 进行 了 有 限 元 模拟 , 得 到 了 压 痕 表面 和 唱 界 
位 置 处 的 应 力 分 布 . 结果 表明 , 相 邻 唱 粒 的 晶体 取 癌 对 压 痕 深度 -载荷 响应 以 及 分 解 切 应 力 分 布 均 有 较 显 著 的 影响 . 压 载荷 
作用 下 , 晶 界 位 置 处 存在 应 力 集中 , 且 唱 体 取向 差异 越 大 , 唱 界 应 力 越 大 . 
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ABSTRACT In the past decades, the indentation test has been widely used to determine the mechanical proper- 
ties of materials. For the micro- or nano- indentation, the indentation response is complex since only one or two 
grains can be indented by the indenter. In order to investigate the indentation behavior of the grain boundary, the 
crystal plasticity theory was implemented into finite element model to simulate the indentation behavior of single 
crystals and bicrystals. The stress distributions on the indented surface and grain boundary were obtained. The re- 
sults showed that the crystallographic orientations of the neighboring grains had a remarkable influence on the 
depth- load response and the resolved shear stress distribution of the indented bicrystals. Under the indentation 
load, stress concentration occurred at the grain boundary, and the stress at the grain boundary increases with the in- 
crease of mis-orientation angle of the two neighboring grains. 
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压 痕 测 试 , 尤其 是 微观 和 纳米 压 痕 , 由 于 其 简 。 压 痕 测试 , 通常 忽略 其 微观 各 向 异性 , 而 在 宏观 上 
单 、 有 效 和 非 破 坏 性 等 优点 , 在 测定 小 体积 材料 力 ”假设 力学 性 能 为 各 向 同性 , 从 而 使 问题 得 到 简化 . 
学 性 能 方面 具有 突出 的 优势 . 对 于 多 蝇 材 料 的 宏观 。 与 此 同时 , 压 痕 尺寸 效应 和 材料 晶体 学 取向 对 纳米 
压 痕 行 为 的 影响 也 得 到 了 广泛 的 研究 呈 . 研究 "" 表 
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RN A ma 明 , 在 纳 观 尺度 内 的 压 痕 响应 与 宏观 有 显著 不 同 . 
收 到 初稿 日 期 :2014-06-25, 收 到 修改 稿 日 期 : 2014-10-13 _. 、 

二 二 微观 和 纳 观 斥 关 斥 涉 通 党 会 斥 存 单个 或 2 个 

作者 简介 , 问 五 柱 , 男 1985 年 生 博 二 对 本 微观 和 纳 玩 里 并 人 人 或 2 了 

DOI: 10.11900/0412.1961.2014.00335 晶 粒 上 . 这 种 情况 下 , 压 痕 行 为 呈现 出 显著 的 各 向 
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异性 特征 , 而 这 些 特征 是 传统 塑性 理论 所 无 法 描述 
的 . 因此 , 基于 晶体 位 错 滑 移 演 化 的 晶体 学 塑性 模 
型 被 广泛 用 于 描述 多 唱 材 料 的 微观 力学 行为 
唱 体 塑性 理论 最 早 可 以 追 亢 到 Taylor7m, Hill 和 
Riceo 的 开创 性 工作 , 随后 晶体 塑性 理论 被 引入 有 限 
元 模拟 中 以 考虑 复杂 边界 条 件 、 塑 性 流动 和 硬化 准 
则 外. 在 过 去 几 十 年 中 , 晶体 塑性 理论 被 成 功 应 用 于 
晶体 材料 力学 行为 的 模拟 叫 忆 . 

在 压 痕 方面 , 晶体 塑性 模型 被 用 来 研究 晶体 取 
癌 对 压 痕 各 向 异性 塑性 行为 的 影响 . 在 单 晶 压 痕 实 
验 研究 方面 , Dong 等 上 9 采用 声 发 射 技术 研究 了 (001) 
和 (110) 取 向 MgO 单 晶 在 压 痕 过 程 中 的 突变 现象 . 
Huang 等 所 研究 了 气相 沉积 一 维 微 晶 须 单 晶 W 纳米 
压 痕 届 服 行为 特性 . 这 些 研究 通过 实验 从 宏观 唯 象 
的 角度 描述 了 单 晶 材料 各 向 异性 的 压 痕 行为 . 另 一 
方面 , 多 唱 材 料 纳 米 压 痕 尺 寸 效 应 和 塑性 各 向 异性 
也 得 到 大 量 研究 . 在 压 痕 尺寸 效应 方面 , 应 变 梯度 
塑性 理论 得 到 广泛 应 用 四. Kang 等 吧 通 过 对 挛 晶 诱 
导 塑 性 (TWIP) 钢 压 痕 实验 研究 了 晶体 取向 对 压 痕 
模 量 和 届 服 的 影响 . 

大 多 数 多 唱 材 料 纳米 或 微观 压 痕 均 为 实验 研 
究 , 难以 获取 纳米 压 痕 细节 信息 , 需要 有 限 元 手段 
作为 辅助 . Xu 等 "9 研究 了 镍 基 单 唱 高 温 合 金晶 体 取 
向 与 压 痕 蠕 变 表面 形 貌 的 相关 性 . Liu 等 中 开展 了 
Cu 单 唱 纳米 压 痕 三 维 唱 体 塑性 有 限 元 模拟 . Li 等 上 m 
开展 了 2024 铝 合金 纳米 压 痕 三 维 晶体 塑性 有 限 元 
模拟 . 然而 先前 的 研究 多 侧重 于 压 痕 响应 , 对 压 痕 
载荷 下 唱 界 的 应 力 分 布 的 研究 鲜 有 报道 . 摸 清 晶 田 
应 力 分 布 规律 是 正确 理解 纳米 压 痕 特 性 的 关键 所 
在 . 本 文 作者 在 前 期 工作 中 中 通过 有 限 元 模拟 比较 
了 不 同 取向 的 fcc 单 唱 合 金 压 痕 蠕 变 和 单 轴 拉 伸 蜂 
变 , 并 建立 了 两 者 之 间 的 联系 , 但 压 痕 载荷 下 品 界 
的 应 力 分 布 尚未 得 到 充分 研究 
本 工作 采用 晶体 塑性 本 构 模 型 对 双 唱 晶 界 压 
痕 开 展 有 限 元 模拟 , 得 到 了 压 痕 面 和 唱 界 应 力 分 布 
特征 . 为 了 便于 对 比 , 对 相应 的 单 晶 压 痕 也 开展 了 
数值 模拟 . 
1 晶体 塑性 本 构 模 型 

本 工作 所 用 模型 为 基于 运动 学 理论 的 晶体 学 
本 构 方程 . 晶体 学 塑性 本 构 方 程 最 先 由 Hill 和 Rice™ 
提出 , 假定 塑性 变形 通过 晶体 位 错 滑 移 和 扩散 变形 
来 实现 , 滑 移 系 分 解 切 应 力 为 滑 移 驱 动力 . 对 于 蝇 
体 学 模型 , 总 变形 梯度 下 可 以 分 解 为 弹性 和 塑性 2 
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下 = 天 (1) 
式 中 ,Fr 表示 唱 格 畸变 和 刚性 转动 所 产生 的 变形 梯 
度 ，F" 则 表示 晶体 沿 着 滑 移 方向 的 均匀 剪 切 所 对 
应 的 变形 梯度 . 
由 于 晶体 变形 和 旋转 , 晶 格 向 量 伸 长 并 旋转 . 
变形 后 滑 移 系 w 上 的 滑 移 方向 可 表示 为 : 
*g _ en 四 O) 
式 中 , m" 为 变形 前 第 w 滑 移 系 滑 移 方向 的 单位 向 
量 , m” 为 晶 格 畸变 后 第 w 滑 移 系 滑 移 方向 的 单位 
向 量 . 晶 格 畸变 后 滑 移 面 法 向 向 量 六 ”可 表示 为 : 
人 = (( FY Yn® (3) 
式 中 , n” 为 变形 前 滑 移 面 的 单位 法 向 量 
当前 状态 的 速度 梯度 工 由 标准 方程 给 出 : 
L=FF "=L+L (4) 
式 中 , 工 为 晶 格 畸变 和 刚性 转动 所 产生 的 速度 梯度 ， 
也 表 示 唱 体 沿 着 滑 移 方向 的 均匀 剪 切 所 对 应 的 速度 
梯度 . 速度 梯度 也 可 表示 为 : 
L=D+W (5) 
式 中 , 刀 和 了 球 分 别 为 延伸 和 自 旋 张 量 . 当 塑 性 变形 
是 由 位 错 滑 移 引起 时 , 塑性 速度 梯度 可 以 表示 为 滑 
移 系 剪 切 率 的 线性 组 合 形式 : 
= > PT (6) 


式 中 ,y" 为 第 a 滑 移 系 滑 移 率 , Schmid 张 量 P* 可 
表示 为 : 


六 (m™ 四 元 娩 十 元 多 四 m™) (7) 


定义 er 为 Cauchy 应 力 张 量 , r 为 加 权 Cauchy 
应 力 张 量 , 则 有 : 


= (8) 

上 (9) 

式 中 , 7 为 滑 移 方 向 上 的 牵引 力 . 分 解 切 应 力 ro 可 
以 通过 Schmid 张 量 表示 为 滑 移 方向 上 的 牵引 力 


分 量 : 


TO = Po .7 (10) 

式 中 ,，P" 为 滑 移 系 a 的 Schmid 因子 . 
对 于 率 相关 的 本 构 方程 , 滑 移 率 可 以 表示 为 分 
解 切 应 力 和 参考 应 力 的 函数 . 滑 移 系 前 切 应 变 率 加 


由 Peirce 等 中 给 出 : 


| 
(Oo (人 ||” 
. (oa) . (a) 
7 加 Yo 网 (1 1) 
8 lg 


式 中 ,，g” 为 参考 前 切 应 力 , m 为 应 变 率 敏感 指数 , 
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7” 为 参考 剪 切 率 . 函数 g" 为 晶体 当前 应 变 硬 化 状 
态 , 假定 a" 仅 与 滑 移 应 变 y 有 关 : 


( 


g =& (Y) (12) 


其 中 ， 
四 
7 


y= > (13) 
式 中 ,y" 为 第 w 滑 移 系 的 滑 移 应 变 . 为 简化 程序 ， 
假定 材料 应 变 硬化 可 由 函数 g" 演化 得 出 


外 py ye (4) 
B 
式 中 , hg 为 y 的 函数 , 在 不 对 p 求 和 时 : 
hgs=qophe (15) 
式 中 , gy 为 潜在 硬化 矩阵 ，j 为 单 硬化 率 ": 
2 Ba 44 
Wh (16) 


式 中 , ho, 7. 和 44 为 滑 移 系 硬化 参数 , 适用 于 所 有 清 
移 系 . 


下 


以 上 本 构 方 程 通过 用 户 材料 子 程序 UMAT 植 
入 ABAQUS 开展 有 限 元 模拟 . 该 程序 的 正确 性 和 有 
效 性 已 经 通过 其 它 相 关 研 究 的 验证 名 
2 模型 建立 

本 工作 所 采用 的 三 维 模型 示意 图 如 图 1 所 示 . 
试 件 为 方形 , 压 头 为 圆柱 形 平 压 头 , 压 头 半径 为 
0.25 mm. 为 了 消除 边界 效应 的 影响 , 压 头 直径 与 试 
件 边 长 之 比 选 为 0.1. 边界 条 件 如 下 : 试 件 底面 约束 = 
向 , 左 侧 约束 向 , 前方 表面 约束 x 方向 (参照 图 1 中 
坐标 系统 ). 考虑 单 晶 (图 1a) 和 双 蝇 (图 1b) 2 种 情 
况 , 在 压 涉 上 方 施加 均 布 力 P=900 MPa. 创建 全 局 华 


况 , 设 定 平面 y=0 为 晶 界 , 将 模型 划分 为 2 个 不 同 取 
向 的 晶 粒 

模型 被 划分 为 6144 个 C3D8R 单 元, 见 图 2. 在 
接近 压 痕 区 域 位 置 处 对 网 格 进行 细 化 . 在 用 户 材料 
子 程序 UMAT 中 , 晶体 取向 采用 Euler 角 控制 . 本 工 
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图 1 压 痕 模 型 示意 图 


Fig.1 Schematic of the indentations performed on single crystal (a) and bicrystals (b) (R 一 radius of the indenter 已 一 


indentation load) 


(a) 


Indented 


图 2 模型 有 限 元 网 格 


Fig.2 Finite element mesh of the model 


(a) mesh of the whole model 


(b) mesh of the indented area 
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作 采 用 的 3 个 基本 晶体 取向 见 表 1. 在 初步 模拟 阶 基 单 晶 合 金 的 压 痕 蠕 变 行为 , 结果 表明 , 压 痕 深度 
段 , 每 个 晶体 取向 被 赋予 整个 模型 , 对 应 了 3 种 单 晶 hh 及 扩展 率 有 具有 以 下 规律 : hioow<hiow<hny, hh woo< 
取向 Case A, Case B 和 Case C. 同时 , 任意 2 个 取向 om< 丈 nn 这 种 现象 一 般 认为 是 由 各 向 异性 所 引起 
的 组 合 产生 3 种 不 同 的 双 唱 取向 Case D, Case E 和 的 不 同 滑 移 系 和 应 力 集中 所 和 致 . 
Case F. 以 上 6 种 情况 及 其 晶体 取向 见 表 1. 一 般 来 说 , 晶体 取向 是 影响 单 晶 力学 性 能 的 如 

本 工作 的 单 晶 和 双 唱 材料 均 采 用 钊 基 单 晶 高 要 因素 , 并 且 取 向 差异 的 大 小 决定 了 不 同 取 问 单 
温 合 金 DD3, 该 材料 被 广泛 应 用 于 燃气 涡轮 叶 睛 力学 性 能 差异 的 大 小 .图 3 给 出 了 3 种 不 同 取向 自 
中 . 表 2 给 出 了 晶 粒 的 力学 属性 , 力学 参数 均 由 钊 基 唱 压 痕 深 度 随 载荷 的 变化 . 值得 注意 的 是 , 压 痕 柔 
单 唱 高 温 合金 DD3 拉 伸 测试 获得 后 . 假定 八 面体 滑 度 (dh/dP) 随 晶体 取向 的 改变 而 变化 , Case A 的 压 痕 


HR 


ml 


I 


了 


移 系 {111}<110> 开 动 , 可 以 模拟 fcc 金 属 在 典型 条 件 ”和 柔 度 最 大 , Case C 的 压 痕 柔 度 最 小 . 由 图 3 可 以 看 
下 的 变形 . 出 , Case A 与 Case C 的 取向 差 角 最 大 , 而 Case B 与 
3 结果 与 分 析 Case C 的 取向 差 角 最 小 . 这 对 于 双 晶 压 痕 响应 的 分 
3.1 晶体 取向 对 讨 痕 深度 扩展 的 影响 析 提 供 了 理论 指导 . 

研究 中 表明, 晶体 取向 对 压 痕 深 度 扩 展 具 有 显 图 4 对 比 了 双 唱 及 相应 组 成 单 晶 的 压 痕 深 度 - 
著 的 影响 . Yan 等 中 研究 了 [001], [011] 和 [111] 取 向 镍 载荷 曲线 . 可 以 看 出 , 双 晶 的 压 痕 和 柔 度 处 于 2 个 组 成 


表 1 晶 粒 Euler 角 
Table 1 Euler angles of the grains 


Crystal Case Grain No. p/() 0/() $b/() 
Single crystal A 1 0.0 0.0 0.0 
B 1 45.0 0.0 0.0 

C 1 65.21 50.76 -63.35 
Bicrystal D 1 45.0 0.0 0.0 
2 0.0 0.0 0.0 

E 1 65.21 50.76 -63.35 
2 0.0 0.0 0.0 

F 1 65.21 50.76 -63.35 
2 45.0 0.0 0.0 


Note: 0 一 precession angle, @—nutation angle, $—rotation angle 


表 2 镍 基 单 晶 高 温 合 金 DD3 力 学 参数 站 


Table 2 Mechanical properties of nickel-based single crystal superalloy DD3™ 


Mechanical parameter Value Unit 
Poisson's ratio v 0.33 
Young's modulus E 210 GPa 
Shear modulus G 83 GPa 
Strain rate sensitivity exponent m 0.02 
Reference shear rate ye 0.03 Ss 
Initial resolved shear stress 1 S00 MPa 
Hardening law parameter ho 600 MPa 
Saturation resolved shear stress 到 585 MPa 
Hardening law parameter AA 1;3 
Reference strain rate Gu 300 MPa 
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单 唱 的 压 痕 柔 度 之 间 , 这 是 由 所 组 成 品 粒 的 取向 和 
力学 性 能 差异 所 引起 . 
2.5 
上 2.0| 要 人 Dp 
去 | 一 上 一 人 3856 NH 
8 > 
3 2 
§ 1.0 上 2 
3 0.5 pa 
a 


200 400 600 800 
P/MPa 


1000 1200 


图 3 晶体 取向 对 深度 -载荷 取向 的 影响 


Fig.3 Effect of crystallographic orientation on the propaga- 


tion of indentation depth with applied load P 
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Fig.4 Depth-load curves of single crystals and bicrystals 
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3.2 分 解 切 应 力 分 布 

在 晶体 学 本 构 模型 中 , 分 解 切 应力 为 晶体 学 滑 
移 的 驱动 力 , 滑 移 系 的 开动 和 请 移 均 与 分 解 切 应 力 
的 大 小 有 关 . 图 $ 给 出 了 不 同 取向 单 晶 压 痕 面 上 的 
分 解 切 应 力 分 布 . 可 以 看 出 , 晶体 取向 对 压 痕 表 面 
分 解 切 应 力 大 小 和 分 布 有 较 大 影响 . 观察 压 痕 附近 
应 力 分 布 特征 可 以 看 出 , 对 于 Case A, Case B 和 
Case C, 应 力 分 布 分 别 近似 呈 四 重 对 称 、 一 重 对 称 和 
三 重 对 称 . 

对 于 纳米 压 痕 , 压 头 有 可 能 会 压 在 唱 界 位 置 
处 . 这 种 情况 下 , 相 邻 唱 粒 属性 差异 对 压 痕 行为 的 
影响 需要 进一步 研究 . 图 6 给 出 了 3 种 不 同 取向 双 
晶 的 压 痕 面 分 解 切 应 力 分 布 . 可 以 看 出 , 在 唱 界 位 
置 处 存在 明显 的 应 力 不 连 续 性 . 相同 压 痕 载荷 下 ， 
Case 卫 所 对 应 的 压 痕 应 力 最 大 , 这 是 由 于 相 邻 2 个 
晶 粒 的 取向 差 角 最 大 . 由 图 6 还 可 以 看 出 , 单个 唱 粒 
的 压 痕 应 力图 案 并 没有 改变 , 仅仅 是 大 小 有 变化 . 
双 唱 压 痕 表面 应 力 分 布 的 形状 是 由 相应 的 单 唱 压 
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图 $ 不 同 取向 单 唱 压 痕 面 上 分 解 切 应 力 分 布 


(a) CasesA,B andD 
(b) Cases A, CandE 
(c) Cases B, C andF 


Fig.S$ Resolved shear stress distributions on the indented 
surface for Case A (a), Case B (b) and Case C (c) 
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痕 分 布 图 拼合 而 成 (图 $), 晶 粒 之 间 的 相互 影响 十 分 图 7 给 出 了 3 个 不 同 取 向 双 电 在 平面 x=0 (垂直 
有 限 . 于 晶 界 ) 上 的 分 解 切 应 力 分 布 . 可 以 看 出 , 在 晶 界 处 

对 于 多 唱 材 料 , 微观 应 力 分 布 与 晶 粒 的 取向 有  ” 存在 明显 的 应 力 不 连 续 现 象 , 同时 在 压 头 下 方 唱 界 
关 . 尤其 是 品 界 处 应 力 分 布 对 相 邻 唱 粒 取向 差 角 较 处 存在 应 力 集 中 . 对 于 Case D 和 Case E, 应 力 相 对 
为 敏感 .Li 等 采用 晶体 塑性 理论 对 2024 铝 合金 开 于 晶 界 呈 非 对 称 分 布 ; 对 于 Case FE, 应 力 相 对 于 唱 界 


展 了 三 维 有 限 元 模拟 , 并 指出 , 较 小 的 唱 界 角 使 得 分 的 . 
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图 6 不 同 取向 双 晶 压 痕 面 上 分 解 切 应 力 分 布 图 7 不 同 取 向 双 晶 平面 x=0 上 的 分 解 切 应 力 分 布 
Fig.6 Resolved shear stress distributions on the indented Fig.7 Distributions of resolved shear stress on the plane of 
surface for Case D (a), Case E (b) and Case F (c) X=0 for Case D (a), Case E (b) and Case F (c) 
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8 不 同 取 向 双 晶 晶 界 0=0) 分 解 切 应 力 分 布 
Fig.8 Resolved shear stress distributions at the grain boundary (=0) for Case D (a), Case E (b) and Case F (c) 
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角 越 小 , 相对 于 晶 界 的 应 力 分 布 越 趋 向 于 对 称 . 
图 8 进一步 给 出 了 压 头 下 方 唱 界 O=0) 上 的 分 
解 切 应 力 分 布 . 可 以 看 出 , 在 坐标 (0.25, 0, -0.06) 和 
(-0.25, 0, -0.06) 处 , 即 压 头 边缘 与 晶 界 的 交点 处 发 
生 严 重 的 应 力 集中 . 相 邻 晶 粒 的 取向 差 角 对 唱 界 处 
应 力 集中 程度 有 显著 影响 . Case 卫 的 唱 界 应 力 最 大 ， 
这 是 由 于 CaseE 所 具有 的 唱 界 角 最 大 . 


4 结论 


(1) 单 晶 压 痕 柔 量 ( 深 度 - 载 答 曲线 斜率 ) 与 晶体 
取向 有 关 . 本 工作 所 研究 的 3 种 取向 的 单 晶 Case A， 
Case B 和 Case C, 其 压 痕 面 上 的 分 解 切 应 力 分 布 分 
别 呈 现 四 重 对 称 、 一 重 对 称 和 三 重 对 称 . 

(2) 相 邻 唱 粒 的 晶体 取向 差 角 会 在 唱 界 位 置 引 
起 应 力 不 连续 和 应 力 集中 , 应 力 分 布 相对 于 唱 界 的 
对 称 性 也 受 晶 体 取向 差 角 影响 . 

(G3) 相 邻 晶 粒 取向 差 角 对 蝇 界 位 置 处 的 应 力 集 
FP 程 度 有 较 大 影响 . 本 工作 中 , Case 王 唱 界 应 力 集中 
最 严重 , 最 大 应 力 位 于 压 头 边缘 与 唱 界 的 交点 处 . 
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